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Abstract. I describe a unique, 20-year-long timing program for the bi- 
nary pulsar B0655+64, the stalwart control experiment for measurements 
of gravitational radiation damping in relativistic neutron-star binaries. 
Observed limits on evolution of the B06554-64 orbit provide new bounds 
on the existence of dipolar gravitational radiation, and hence on violation 
of the Strong Equivalence Principle. 


1. Introduction 

PSR B0655+64, in a highly circular one-day orbit with a ~ 0.8 M@ white-dwarf 
companion, serves as a control experiment for measurements of orbital decay 
in the highly relativistic double-neutron-star binaries: General Relativity (GR) 
has predicted equally well the strong back-reaction to gravitational radiation for 
PSRs B1913+16 (Taylor 1993; see also Weisberg & Taylor, this volume) and 
B 1534+12 (Stairs et al. 1998), and the absence of detectable orbital evolution 
for PSR B0655+64 over two decades. The long-term stability of the B0655+64 
orbit sets unique bounds on departures from GR that give rise to dipolar grav- 
itational radiation (Arzoumanian 1995, Goldman 1992), the existence of which 
would represent a violation of the Strong Equivalence Principle (SEP), one of 
the basic tenets of GR. Gerard & Wiaux (2002) examine the theoretical basis 
for dipolar gravitational radiation and suggest that bounds from binary pulsars 
may be competitive with future satellite experiments dedicated to probing SEP 
violation. 

A definitive analysis of the available data for PSR B0655+64, and impli- 
cations for a variety of alternative theories of gravitation, will be published 
separately (Arzoumanian et al. 2003, in preparation). Following recent observa- 
tions to extend our long-term monitoring program, I present here preliminary 
results on the orbital evolution of PSR B0655+64. 

1.1. Binary Pulsars and the Strong Equivalence Principle 

The Strong Equivalence Principle posits that the response of a body to an ex- 
ternal gravitational field is independent of the gravitational self-energy of the 
body; it would be violated if, for example, two objects with different gravita- 
tional binding energies were observed to fall at different rates. SEP is satisfied 
by postulate within GR; as a consequence, the lowest allowed multipole order 
of gravitational radiation is the “electric” quadrupole. If SEP does not hold, 
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however, emission of dipolar gravitational radiation is allowed, and lends itself 
in principle to indirect detection through decay of a binary orbit at a rate incon- 
sistent with the quadrupole-order prediction of GR. Because self-gravitational 
effects are important in neutron stars, SEP can be tested with pulsar systems. 

The most thoroughly studied theory of gravitation that violates SEP, the 
so-called “Brans-Dicke” scalar-tensor theory (Brans & Dicke 1961), invokes the 
existence of a scalar gravitational field that couples to matter through a single di- 
mensionless parameter, cvbd- The Brans-Dicke theory becomes indistinguishable 
from GR as wbd becomes arbitrarily large; solar-system experiments currently 
place a lower bound u>bd ^ 500. Other alternatives to GR are also constrained 
by limits on dipolar gravitational radiation (see, e.g., Damour & Esposito-Farese 
1992): recent observations of distant type la supernovae suggesting that the ex- 
pansion of the universe is accelerating have prompted renewed interest in theories 
involving scalar fields (e.g., “quintessence,” Caldwell et al. 1998) or long-range 
forces that may be responsible for the acceleration. Such fields would radiate 
predominantly at dipole order. 

The various contributions to observed changes in the orbital period Pb of a 
binary system (neglecting tidal and mass-transfer effects) can be expressed as 



where the subscripts Q, D, G, and a denote the effects of quadrupolar and 
dipolar gravitational radiation, a change in the gravitational constant with time, 
and the varying Doppler shift from a relative acceleration of the solar system and 
pulsar binary 1 . In a large class of metric theories of gravitation, the radiation 
terms for a circular orbit are given by (Eardley 1975, Will 1981, Goldman 1992) 


jQ ~ 4/ 3 (§) ( 2 ) 

k d Q~ 2 (si - s 2 ) 2 n 2 ^p-T Q , (3) 

where to = 2n/Pb, mi and m 2 are the pulsar and companion masses in solar units, 
M = mi + 7712 , and T© is the mass of the Sun in units of time. The dimensionless 
parameters «i, up, and Q describe the strength of quadrupolar and dipolar 
radiation and the effective gravitational constant in a given theory. In GR, 
Ki — 12, kd = 0, and Q = 1; in the Brans-Dicke theory, kdQ~ 2 = 2(2 + wbd) - 1 . 
At the level of the current constraint on the dipole term D, the G contribution 
to Pb in Eq. 1 can be neglected (Arzoumanian 1995, Kaspi et al. 1994). 

The quantities Si and S 2 in Eq. 3 represent the “sensitivities” of the orbiting 
objects, the fractional change in binding energy of each star with G , 




'To correct for Galactic accelerations and the “Shklovskii” effect, I follow Damour & Taylor 
(1991), using a new proper-motion measurement for B0655+64: n = 6.8 ± 1.1 masyr -1 . 
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where N is the total number of baryons. The sensitivity of white dwarfs is negli- 
gible; neutron stars are thought to have sensitivities in the range 0.15 < s < 0.40 
(Will & Zaglauer 1989), with larger values corresponding to “softer” equations 
of state. Results of X-ray burst oscillation modeling (e.g., Nath et al. 2002), 
and the recent detection of atmospheric absorption lines from a neutron star 
(Cottam et al. 2002), suggest that fairly stiff equations of state are appropriate. 
I therefore adopt a fiducial value s = 0.2 below. It is worth noting that the 
difference (si — S 2) 2 in Eq. 3 strongly suppresses dipolar radiation from NS-NS 
binaries like B1913+16, so that the celebrated agreement of the latter’s orbital 
decay rate with the GR prediction does not usefully constrain the existence of 
gravitational radiation at dipole order. In scalar-tensor theories, gravitational 
binding energy takes on the role of an effective gravitational “charge.” Neu- 
tron star-white dwarf (NS-WD) binaries can therefore be expected to copiously 
emit dipolar radiation, if such radiation exists, because of the large difference in 
gravitational self-energy between the component stars. 

1.2. PSR B0655+64 

PSR B0655+64 was discovered during a survey of the Northern sky made with 
the NRAO 300 Foot telescope (Damashek et al. 1982). Timing observations 
began soon thereafter, and the companion was optically identified as a massive 
white dwarf by Kulkarni (1986). If we assume pulsar and companion masses of 
1.4 Mq and 0.8 M@, Eq. 2 predicts a rate of orbital decay within GR, P GR = 
—2 x 10 -14 , corresponding to a decay timescale for the orbit of tq ~ 150 Gyr. 


2. Observations and Results 

The current data span 20 years, including observations made with the NRAO 
Green Bank 300 Foot (Taylor & Dewey 1988), 140 Foot (Backus et al. 1982, 
Arzoumanian et al. 1994b), and GBT telescopes, as well as the NRAO VLA 
(Thorsett 1991) and Jodrell Bank Lovell (Jones & Lyne 1988) telescopes. We 
carried out intensive observing campaigns, to obtain good coverage of all orbital 
phases at a single epoch, every few years with the 140 Foot and recently with the 
GBT, and these data provide the interesting constraints on orbital evolution. 

Pulse times-of-arrival (TOAs) were derived from observations following stan- 
dard procedures and analyzed with the Tempo software package, using the ELLl 
binary timing formula (Lange et al. 2001). A least-squares fit of the entire 
dataset for rotational, astrometric, and orbital parameters produces residuals 
consistent with statistical fluctuations. 

Orbital phase residual and period measurements over time are shown in 
Figure 1. The constancy of Pj, confirms the result of the global timing solu- 
tion: the available data bound changes in orbital period at the level |Pft| 0 bs < 
1.5 x 10 -13 (la), a limit ~ 7 times the GR prediction. Accounting for a small 
correction for relative acceleration (Eq. 1), we then have a 2a upper limit on 
orbital evolution, 

< 1.0 x KT 10 yr -1 (2a). 
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Figure 1. Orbital evolution of PSR B0655+64. Box outlines depict 
ranges of dates over which phase and period measurements were made. 


The resulting constraints on u>bd are summarized in Fig. 2. Horizontal 
lines depict the current la and 2 a limits on orbital decay in the PSR B0655+64 
system. The curving dashed and dotted lines are the expected orbital period 
changes due to emission of gravitational radiation through quadrupole order, 
Eqs. 2-3, for NS sensitivities from soft and stiff equations of state respectively. 
The three curves in each set represent the (mi, m 2 ) pairs (1.30, 0.7), (1.35, 
0.8), and (1.40, 0.9), so that total system mass increases to the upper left. 
Constraints on ubd then lie at the intersections of the dashed or dotted curves 
with the measured upper limits on P&. From Eqs. 2, 3, and 5, we have in general 

k d G- 2 < 0.006 (^) 2 (2a), (6) 

and for the Brans-Dicke theory specifically, ojbd it 320 (s/0.2) 2 (2a). 


3. Discussion 

Prospects for improving bounds on SEP violation and the existence of dipolar 
gravitational radiation are good. Damour & Taylor (1992) show that, for con- 
stant data quality, measurement uncertainty for Pt> scales with data span T as 
T -5 / 2 ; moreover, data quality typically improves with time through improved 
instrumentation. Also, pulsar surveys continue to discover NS-WD systems sim- 
ilar to B0655+64 but with shorter orbital periods and millisecond pulse periods. 
Their higher rates of gravitational energy release coupled with higher timing 
precision suggest that these systems will surpass B0655+64 as laboratories for 
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Figure 2. Bounds on ubd from PSR B0655+64 (after Will & Za- 
glauer 1989). 


testing theories of gravitation — Table 1 lists the current sample of relativistic 
NS-WD systems. While B0655+64 remains a valuable object of further study, a 
potential caveat applies to relativistic studies of all close NS-WD binaries: very 
small orbital separations can give rise to non-relativistic interactions. Heating of 
the companion star by irradiation from the pulsar is thought to power tides and 
mass loss in a handful of low-mass binaries. Outflows and tides cause significant 
orbital torques (e.g., Arzoumanian et al. 1994a), which would overwhelm any 
small secular trend due to gravitational radiation. If indeed the closest NS-WD 
binaries (e.g., PSR B075H-18; see Nice et al., this volume) are “clean” systems, 
significant new constraints on SETP violation will emerge in the coming years. 
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Table 1. Relativistic NS-WD systems. Spin and orbital periods are 
shown alongside spin-down luminosity (E), orbital separation (/l) and 
timescales for orbital evolution from quadrupolar ( tq ) and dipolar (to) 
radiation. The potential for heating and tidal interactions between the 
component stars increases rapidly for increasing E and decreasing A. 
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